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Abstract. Three-dimensional magnetohydrodynamic simulations of the sur- 
face layers of the Sun intrinsically produce a predominantly horizontal magnetic 
field in the photosphere. This is a robust result in the sense that it arises 
from simulations with largely different initial and boundary conditions for the 
magnetic field. While the disk-center synthetic circular and linear polarization 
signals agree with measurements from Hinode, their center-to-limb variation sen- 
sitively depends on the height variation of the horizontal and the vertical field 
component and they seem to be at variance with the observed behavior. 



1. Introduction 



Observations with the spectropolarimeter of the Solar Optical Telescope (SOT) 
onboard the Hinode space observatory (Kosugi et al. 20071 indicate that seen 



with a spatial resolution of 0.3", quiet internetwork regions harbor a photo- 
spheric magnetic field whose mean field strength of its horizontal component 
considerably surpa sses that of the vertical component ( Lites et al.| 2008 Orozco 
Suarez et al. 2007 1 . According to these papers, the vertical fields are concen- 



trated in the intergranular lanes, whereas the stronger, horizontal fields occur 
most commonly at the edges of the bright granules, aside from the vertical fields. 
Harvey et al. (2007) find from recordings with GONG and SOLIS a 'seething 
magnetic held' with a line-of-sight component increasing from disk center to 
limb as expected for a nearly horizontal field orientation. Ishikawa et al. (20081 



detected transient horizontal magnetic fields in plage regions as well. Previously, 
Martinez Fillet et al. (19971 and Meunier et al. (1998) reported observations of 



weak and strong horizontal fields in quiet Sun regions. 



Regarding numerical simulations, Grossmann-Doerth et al. (1998), note "we 



find in all simulations also strong horizontal fields above convective upflows" , and 
Schaffenberger et al. (2005, 2006) find frequent horizontal fields in their three- 
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Figure 1. Horizontally and temporally averaged horizontal (solid) and ver- 
tical (dashed) absolute field strength as functions of height for run h20 (heavy 
curves) and for run vlO (light curves). 



dimensional simulations, which they describe as "small-scale canopies". Also 
the 3-D simulations of Abbett (20071 display "horizontally directed ribbons of 
magnetic flux that permeate the model chromosphere", not unlike the figures 
shown by ^Schaffenberger et al. ( 2006 1 . [Schiissler &: Vogler ( 2008| ) find in a three- 
dimensional surface-dynamo simulation "a clear dominance of the horizontal 
field in the height range where the spectral lines used for the Hinode observations 
are formed" . 



Steiner et al. (20081 report on results from three-dimensional magnetohy- 



drodynamic numerical simulations of the internetwork magnetic field with regard 
to the intrinsically produced horizontal magnetic field. They compute the po- 
larimetric signal of this field and compare it to measurements with Hinode. In 
the following, we briefly summarize part of their results and present new cal- 
culations of the center-to-limb variation of the linear and circular polarization 
from these simulations. 



2. Simulations 



We have carried out three sets of numerical simulations, run vlO, run h20, and 
run h50, which significantly differ in their initial and boundary conditions for 
the magnetic field. Run vlO starts with a homogeneous vertical field of 10 G 
strength, which is kept vertical at the top and bottom boundary. These condi- 
tions might actually be more appropriate for network magnetic fields because of 
the preference for one polarity and the vertical direction. Run h20 and run h50 
start without a magnetic field but upwellings that enter the simulation domain 
across the bottom boundary area carry horizontal magnetic field of a uniform 
strength of 20 G and 50 G, respectively, into the box. The three-dimensional 
computational domain extends from 1400 km below the mean surface of optical 
depth Tc = 1 to 1400 km above it, high enough that the top boundary condition 
should not unduly tamper the formation layers of the spectral lines that are used 
in polarimetric measurements with Hinode. 



The horizontal magnetic Geld 
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Figure 2. Center-to-limb variation (right to left) of the total (wave- 
length integrated) circular (T4ot) and linear (itot) polarization in pni of the 
Fe I 630.25 nm spectral line of a snapshot of run vlO (left) and h50 (right). 



Figure [T] shows the horizontally and temporally averaged absolute vertical 
and horizontal magnetic field strength as functions of height of the magnetic 
field, which emerges from the two runs vlO and h20. In run h20, the mean hori- 
zontal field strength is larger than the mean strength of the vertical component 
throughout the photosphere and the lower chromosphere: in run vlO this is the 
case in the height range between 250 km and 850 km. Both runs show a local 
maximum in the horizontal field strength close to the classical temperature min- 
imum at a height of around 500 km. At a given continuum optical depth in the 
photosphere the area occupied by fields with a horizontal component stronger 
than a given threshold of a few gauss is typically three times larger than the 
area with a vertical component exceeding this limit. 

For a comparison with Zeeman measurements from the Hinode spectropo- 
larimeter we synthesized the Stokes profiles of both 630 nm Fe i spectral lines 
of the two simulation runs. Profiles were computed with the radiative transfer 
de SIR ( Ruiz Cobo Sz del Toro Iniesta|1992 1 with a spectral samplin g of 2 pm. 



CO 



When applying an appropriate point spread function, PSF, (Wedemeyer-Bohm 



20081 to the synthetic profiles and subjecting them to the same procedure for 



conversion to apparent fiux densities as done by Lites et al. (20081 for the ob- 
served profiles, we obtain spatial and temporal averages for the transversal and 
longitudinal apparent magnetic fiux densities, -Bjpp and |-Bappl of respectively 
21.5 Mxcm~^ and 5.0 Mxcm~^ for run h20 and 10.4 Mxcm"^ and 6.6 Mxcm~^ 
for run vlO. Thus, the ratio r = (5jpp)/(|Sappl> = 4.3 for h20 and 1.6 for 



Lites et al. 



vlO. 



(I^appl) = 11 Mxcm 



(20081 obtain from Hinode SP data 
resulting in r 



app 



I) = 55 Mxcm and 



5.0. 



3. Center-to- limb variation (CLV) 

Figure [2] shows the center-to-limb variation of the wavelength integrated absolute 
circular (Vtot) and absolute linear (Ltot) polarization signal, for the 630.25 nm 
Fe I spectral line in pm. The left panel of Fig. |2]shows the result for run vlO, the 
right panel for run h50, both for one single simulation snapshot. For these results 
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foreshortening was taken into account but no PSF. The computational domain 
of h50 is in the horizontal directions twice as large as that for vlO and h20. The 
CLV of Vtot and Ltot substantially differ between vlO and h50, which indicates 
that they sensitively depend on the precise height variation of the horizontal 
and the vertical component of the magnetic field. Since the vertical component 
of vlO dominates over the horizontal one in the low photosphere (cf. Fig. [T]), the 
mean longitudinal (line-of-sight) field component decreases when moving away 
from disk center and only when approaching the limb it increases again due to 
the sharp increase of the horizontal component above z ~ 200 km. This behav- 
ior is reflected in the CLV of Vtot of run vlO. The transversal component, viz. 
the component perpendicular to the line-of-sight and hence Ltoti increases at 
first because the vertical component stays almost constant with height in model 
vlO (in contrast to model h20 or h50). For h50 we have an increase in the longi- 
tudinal component as a function of disk-center distance because the horizontal 
component dominates throughout the photosphere and as a consequence, the 
transversal component decreases. This behavior is reflected by the CLV of Vtot 
and Ltot of run h50. 

While the CLV of Fig. |2] finds a natural explanation in terms of the ver- 
tical stracture of (-Bhor) and (-Bvcr) shown in Fig. [T] the observations by Lites 



et al. (2008) tell a different story. Accordingly, the circular polarization signal 
decreases from the disk center to the limb as it does in Fig. [2] (left). But also 
the linear polarization signal decreases like in Fig. [2] (right). None of the two 
simulation snapshots to which Fig. [2] refers show a simultaneous and monotonic 
decrease in both quantities. 
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